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Abstract: Gene expression changes in response to aging, heat stress, hyperoxia, hydrogen peroxide, and ionizing radiation 
were compared using microarrays. A set of 18 genes were up-regulated across all conditions, indicating a general stress 
response shared with aging, including the heat shock protein (Hsp) genes Hsp70, Hsp83 and l(2)efl, the glutathione-S- 
transferase gene GstD2, and the mitochondrial unfolded protein response (mUPR) gene ref(2)P. Selected gene expression 
changes were confirmed using quantitative PCR, Northern analysis and GstD-GFP reporter constructs. Certain genes were 
altered in only a subset of the conditions, for example, up-regulation of numerous developmental pathway and signaling 
genes in response to hydrogen peroxide. While aging shared features with each stress, aging was more similar to the 
stresses most associated with oxidative stress (hyperoxia, hydrogen peroxide, ionizing radiation) than to heat stress. Aging 
is associated with down-regulation of numerous mitochondrial genes, including electron-transport-chain (ETC) genes and 
mitochondrial metabolism genes, and a sub-set of these changes was also observed upon hydrogen peroxide stress and 
ionizing radiation stress. Aging shared the largest number of gene expression changes with hyperoxia. The extensive 
down-regulation of mitochondrial and ETC genes during aging is consistent with an aging-associated failure in 
mitochondrial maintenance, which may underlie the oxidative stress-like and proteotoxic stress-like responses observed 
during aging. 



INTRODUCTION 



Heat shock protein (Hsp) genes are induced in response 
to stresses that cause protein denaturation, through 
activation of the heat shock factor (HSF) [1]. Up- 
regulation of Hsp genes is also observed during normal 
aging [2], For example, both Hsp70 and Hsp22 are up- 
regulated during normal Drosophila aging, and this up- 
regulation requires functional HSF binding sites (Heat 
Shock Elements, or HSEs) in the promoters of these 
genes [3-5]. Genome-wide studies of gene expression 
changes during Drosophila aging have revealed 
additional features of a stress response, including the 
up-regulation of additional oxidative stress-response 



genes, and the dramatic up-regulation of innate immune 
response genes [6-8]. In addition, Drosophila aging is 
characterized by a small but across-the-board down- 
regulation of mitochondrial metabolism and electron 
transport chain (ETC) genes [6, 8], and this pattern is 
also observed in aging mammalian tissues [9], and at 
early adult ages in both Drosophila and C. elegans [10], 
indicating a conservation of aging mechanisms across 
species. Both innate immune response genes [6] and 
Hsp genes [11, 12] have been shown to be predictive 
biomarkers of individual animal life span when the gene 
promoters are fused to GFP to create transgenic 
reporters, thereby supporting the significance of the 
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identified gene expression changes. Here normal aging 
was compared with multiple stressors to provide further 
insight into common and unique features. 

RESULTS 

Gene expression changes common to each stress and 
to aging 

Micro-array analysis was used to identify genes whose 
expression was altered in response to normal aging, 
hyperoxia, hydrogen peroxide, ionizing radiation and 
heat stress. A core set of 18 stress-response genes were 
up-regulated >1. 5-fold in response to each of the tested 
stresses as well as during normal aging (Table 1). 



These up-regulated genes included the heat shock 
protein genes Hsp70, Hsp83 (which is the single 
Drosophila Hsp90-class member), and the small Hsp 
gene l(2)efl. The up-regulation of Hsp70 and 1(2) efl in 
response to selected stressors was confirmed using 
quantitative real-time PCR analysis (Figure 1), and in 
addition Hsp70 was analyzed by Northern blot analysis 
(Supplemental Figure SI; results summarized in Table 
2). Also up-regulated by aging and each stressor were 
the glutathione S-transferase gene GstD2, the central 
metabolic regulatory enzyme gene Pepck, and the 
mitochondrial unfolded protein response (mUPR) gene 
ref(2)P. Down-regulated genes included several 
associated with sugar metabolism and proteolysis 
(Table 1). 
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Figure 1. Hsp70 and l(2)efl RNA levels in response to selected stresses. Quantitative real-time RT-PCR analysis 
was used to determine RNA levels for the genes Hsp70 and l(2)efl in response to selected stresses, in both male and 
female flies, as indicated. HS, heat stress; IR, ionizing radiation. Stress treatment RNA levels were compared to control 
using unpaired, two-sided t-tests, and statistically significant differences (p < 0.05) are indicated with asterisk. 
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Table 1. Gene expression changes common to aging and 
each stress 

a) 18 genes up-regulated in aging and all other stresses 



CG6489 


Hsp70 


Heat-shock-protein-70 


CG3705 


aay 


astray 


CG32130 


stv 


starvin 


CG4533 


l(2)efl 


lethal (2) essential for life 


CG33229 


CG33229 




CG4181 


GstD2 


Glutathione S transferase D2 


CG3821 


Aats-asp 


Aspartyl-tRNA synthetase 


CG5966 


CG5966 




CGI 1030 


CGI 1030 




CGI 4245 


CGI 4245 




CG14246 


CG14246 




CG15784 


CG15784 
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CG31638 




CG13941 


Arc2 


Arc2 


CG1242 


Hsp83 


Heat shock protein 83 


CG10360 


ref(2)P 


refractory to sigma P 


CG32103 


CG32103 




CG17725 


Pepck 


Phosphoenolpyruvate carboxykinase 



b) GO enrichment terms for genes upregulated in aging and all other stresses 
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heat shock-mediated polytene chromosome puffmg(5) 
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response to heat(7) 
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response to temperature stimulus(7) 


5.70E-07 
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cellular response to heat(5) 


1.12E-06 
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response to oxygen levels(5) 
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2.00E-04 



www.impactaging.com 



770 



AGING, November 2012, Vol.4 No.ll 



c) 32 genes down-regulated in aging and all other stresses 



CG number 


Symbol 


Gene name 


CG10026 


CG10026 
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CGI 0467 
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CG14120 




CGI 4661 


CG14661 
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CG18302 




CGI 8493 


CGI 8493 
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CG5107 
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CGI 1669 


Mal-A7 


Maltase A7 
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d) GO enrichment terms for genes down-regulated in aging and all other str esses 
GO:0006013 mannose metabolic process(3) 0.018117 
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Table 2. Confirmation of selected gene expression changes using 
qPCR and Northern analysis 
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XQL 


XN 
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XQL 
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XQ 


XQ 


XL 


Drs 
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X Q 


XQf 


XL 


ade3 


XQL 


Xnt 


Q 


XQ 


XL 


CGI 1089 


XQ 


XN 


Q 


XQ 


Xnt 




X GeneChip data (this study); Q q-PCR analysis (this study); N Northern 
analysis (this study); L Northern analysis (Landis et al 2004 PNAS 101:7663-8); 
nt not tested. a fold increase >1.2 



Gene expression changes unique to each stress 

Each stress had gene expression changes that were 
unique to that stress (listed in Supplemental Table SI) 
and the enriched GO terms that uniquely characterize 
each stress are summarized (Table 3). Hyperoxia stress 
had no enriched GO terms in the uniquely up-regulated 
genes, and a single enriched GO term, Signal peptide 
processing (3 genes) among the down-regulated genes. 
In contrast, there were numerous up-regulated genes 
unique to hydrogen peroxide stress, and these up- 
regulated genes were enriched for many GO terms 
involved in developmental pathways, signaling 
pathways, and nucleobase metabolism (Table 3). Genes 
uniquely up-regulated upon heat stress included many 
of the Hsp60-class, and this list was consequently 
enriched for the GO term Protein folding (16 genes), 
whereas down-regulated genes unique to heat stress 
were enriched for the GO terms Defense response and 
Melanization defense response (Table 3). Finally, 
genes uniquely up-regulated in response to ionizing 
radiation included several proteasome subunit genes 
(Supplemental Table SI), and this gene list was 
enriched for the GO terms Protein catabolic process and 
Macromolecular catabolic process (Table 3), whereas 
there were no GO terms enriched among down- 
regulated genes. 

Aging is most similar to hyperoxia 

As described above, a core set of stress response genes 
was induced during aging and by each of the stressors 
tested. Aging shared additional changes in gene expres- 



sion with each individual stressor (Supplemental Table 
S2), and was found to be more similar to the stresses 
most associated with oxidative stress (hyperoxia, 
hydrogen peroxide, ionizing radiation) than it was to 
heat stress, based on cluster analysis (Supplemental 
Figure S2) and by comparison of the GO categories that 
were enriched in the groups of up-regulated and down- 
regulated genes (Supplemental Table S3). While aging 
shared a significant overlap in up-regulated and down- 
regulated genes with each of the stresses, aging shared 
the greatest number of gene expression changes with 
hyperoxia (Table 4). 

Gene expression changes unique to aging 

A number of gene expression changes were found to be 
unique to aging. These included up-regulation of 
numerous innate immune response genes, and down- 
regulation of numerous mitochondrial metabolism 
genes, including ones encoding components of the ETC 
(Supplemental Table SI; enriched GO terms listed in 
Table 4). While up-regulation of innate immune 
response genes is a feature of aging that is shared with 
hyperoxia [6] (Supplemental Table S3), the number of 
up-regulated innate immune response genes was 
significantly greater for aging, resulting in many 
changes in this category that were unique to aging. 
Also uniquely up-regulated during aging were the 
odorant receptor genes Obp56a and Obp57d. 

Down-regulation of mitochondrial genes is a feature of 
aging that is shared with hydrogen peroxide and 
ionizing radiation (Supplemental Table S3), but the 
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number of down-regulated mitochondrial genes was mitochondrial ribosomal proteins and components of the 

greater for aging, resulting in many changes in this mitochondrial membrane protein translocases (TIM and 

category that were unique to aging. Among these many TOM), as well as the mitochondrial form of superoxide 

down-regulated mitochondrial genes were ones encoding dismutase (MnSOD or Sod2). 



Table 3. Features unique to each stress 

a) GO enrichment terms for genes uniquely up-regulated in hyperoxia 
None found. 



b) GO Enrichment Terms for genes uniquely down-regulated in hyperoxia 
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signal peptide processing(3) 
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GO:0050896 


response to stimulus(106) 


0.001425 


GO:0051239 


regulation of multicellular organismal process(34) 
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GO:0045595 


regulation of cell differentiation^ 5) 


0.001776 


GO:0030182 


neuron differentiation(39) 


0.00189 


GO:0023052 


Signaling(79) 


0.002054 


GO:0007154 


cell communication(80) 


0.002305 


GO:0048666 


neuron development(35) 


0.003766 


GO:0007165 


signal transduction(64) 


0.00438 


GO:0048513 


organ development(62) 


0.004542 


GO:0022414 


reproductive process(54) 


0.005853 


GO:0003002 


Regionalization(3 3 ) 


0.007295 


GO:0022603 


regulation of anatomical structure morphogenesis(21) 


0.00904 


GO:2000026 


regulation of multicellular organismal development^ 6) 


0.01046 


GO:0009880 


embryonic pattern specification^ 1) 


0.012338 


GO:0009887 


organ morphogenesis(37) 


0.013194 



www.impactaging.com 



773 



AGING, November 2012, Vol.4 No.ll 



00:0007350 


blastoderm segmentation(20) 


0.01346 


GO:0003006 


developmental process involved m reproduction(37) 


A A 1 O A A O 

0.018008 


GO:0051093 


negative regulation of developmental process(16) 


A A 1 O O A 1 

0.018841 
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regulation of macromolecule biosynthetic process(49) 


A A 1 A 1 AO 
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A A 1 A 1 no 
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d) GO enrichment terms for genes uniquely down-regulated in hydrogen peroxide 




None found. 






e) GO enrichment terms for genes uniquely upr-regulated in heat stress 




GO:0006457 


protein folding(16) 


0.018356 


f) GO enrichment terms for genes uniquely down-regulated in heat stress 
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melanin metabolic process(7) 


0.001782 
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melanization defense response(6) 
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GO:0006952 
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g) GO Enrichment terms for genes uniquely up-regulated in Ionizing radiation 
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h) GO enrichment terms for genes uniquely down-regulated in Ionizing radiation 
None found. 
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Table 4. Number of gene expression changes shared by aging and individual stresses 



Aging change 


Number genes 


Stress change 


Number genes 


Number in 
common 
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Aging up 


456 


Hyperoxia up 


335 


165 


9.8 x 10 -170 


Aging up 


A ^£ 


Ionizing 
radiation up 


/ 10 


i/i 


A.y X 1U 


Aging up 


A ^£ 


Hydrogen 
peroxide up 


/ Zo 


1 ii 
i j j 


a i v 1 fr^o 

4. J X 1U 


Aging up 


A^& 


Heat stress up 


/ J4 


DO 


1 1 v 1 PT 15 
l.ZX 1U 


Aging down 


1009 


Hyperoxia 
down 


556 


222 


1.5 xlO" 121 


Aging down 


1009 


Ionizing 
radiation down 


674 


166 


3.6 x 10" 54 


Aging down 


1009 


Hydrogen 
peroxide down 


911 


132 


1.1 x 10" ls 


Aging down 


1009 


Heast stress 
down 


806 


89 


7.3 x 10" 7 



Confirmation of selected gene expression changes 

Changes in gene expression caused by one or more 
stressors were confirmed by quantitative real-time PCR 
(Figure 1 and Supplemental Figure S3) and by Northern 
blot analysis (Supplemental Figure SI), and in general 
an excellent concordance was observed with the micro- 
array data and with the published literature 
(Summarized in Table 2). One exception was for 
expression of the innate immune response gene 
Drosomycin upon heat stress, which was observed to 
increase in the micro-array analysis, but to decrease in 
the qPCR analysis (Table 2). Because bacterial load and 
Drosomycin gene expression can vary significantly 
between different flies and vials of flies [13], we 
conclude that this discrepancy was most likely due to a 
small difference in bacterial load and Drosomycin gene 
expression in the control flies used for the qPCR 
analysis relative to the control flies used for micro-array 
analysis. 

Comparison of the responses to the different stresses 
reveals preferential induction of certain genes. For 
example, Hsp70 (Figure 1) and Hsp22 (Supplemental 
Figure S3) were induced to the greatest extent by heat 
stress, whereas l(2)efl (Figure 1) and ade3 
(Supplemental Figure S3) were induced to a greater 
extent by ionizing radiation and hyperoxia. In addition 
significant sexual dimorphism in the magnitude of 
responses was observed. For example, the induction of 
Hsp70 (Figure 1) and Hsp22 (Supplemental Figure S3) 
in response to heat stress was greater in males than in 
females, and the induction of 1(2) efl (Figure 1) and ade3 



(Supplemental Figure S3) in response to ionizing 
radiation was greater in males than in females. 

A GstD-GFP reporter construct recapitulates 
induction during aging 

The GstDl gene encodes a glutathione-S-transferase, 
and is induced in adult flies during normal aging and 
when flies are challenged with oxidative stress 
produced by hyperoxia and paraquat [6, 7, 14], and was 
also found to be up-regulated in response to hydrogen 
peroxide and ionizing radiation stress (Summarized in 
Table 2). The GstDl promoter region contains 
consensus binding motifs for the stress-responsive 
transcription factors Nrf2 and Foxo (diagrammed in 
Supplemental Figure S4). A transgenic reporter has 
been characterized where the regulatory sequences of 
the GstDl gene are fused to GFP, and the resulting 
GstD-GFP reporter is induced in the adult fly by 
feeding flies with the oxidative stressors paraquat, 
arsenic or hydrogen peroxide [15]. A clustered point 
mutation was created to disrupt the antioxidant response 
element (ARE) in the GstD-GFP reporter to yield a 
mutant reporter called GstD-deltaARE-GFP 
(diagrammed in Supplemental Figure S4). These 
reporters have been used to demonstrate that the GstD- 
GFP transgene is positively regulated in the adult fly in 
response to genetically-altered Nrf2 expression, and in 
response to the cancer chemotherapeutic compound 
Oltipraz which is known to activate Nrf2, and these 
regulations required the intact ARE [15]. Quantitative 
PCR analysis of adult flies indicated that induction of 
the GstDl gene by paraquat is reduced in flies hemizy- 
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gous for the JNKK gene hemipterous, suggesting 
additional positive regulation of GstD genes by the JNK 
pathway in response to oxidative stress [16]. The JNK 
pathway activates the transcription factor Foxo 
suggesting that the JNK pathway may activate GstD 
gene expression through the Foxo binding motif located 
in the GstD gene promoter region [17] (diagrammed in 
Supplemental Figure S4). A GstDl-LacZ reporter has 
been reported to be up-regulated during normal aging in 
the enteroendocrine cells (ECs) of the fly intestine, but 
to be reduced during aging in the intestinal stem cells 
(ISCs)[18]. 

GstD-GFP GstD-deltaARE-GFP 




Figure 2. GstD-GFP transgenic reporters recapitulate GstDl 
gene induction during aging. (A) Expression of the transgenic 
reporter constructs GstD-GFP and GstD-deltaARE-GFP was visualized in 6 
day old (Young) and 30 day old (Old) male and female flies, as indicated, 
using the fluorescence stereomicroscope. The GFP image and an 
overlay of the GFP image and the visible light image are presented, as 
indicated. (B) Quantification of the expression of the GstD-GFP and 
GstD-deltaARE-GFP reporters in male (M) and female (F) flies, both 
young (Y) and old (0), as indicated. The data for the Gst-deltaARE-GFP 
reporter is specified by prefix (delta). Values for old flies were 
compared to young using unpaired, two-sided t-tests, and statistically 
significant differences (p < 0.05) are indicated with asterisk. 



Here the expression of the GstD-GFP and GstD- 
deltaARE-GFP reporters were examined in whole adult 
flies during normal aging. The GstD-GFP reporter was 
expressed at low levels in young flies, and exhibited 
robust induction throughout the body of the fly during 
aging, including the head, thorax, abdomen and legs 
(Figure 2A), consistent with the whole-body micro- 
array analyses presented above. Induction was apparent 
even at moderate ages (30 days; Figure 2A) and 
continued at high levels for the remainder of the life 
span (data not shown). In contrast, the GstD-deltaARE- 
GFP reporter was robustly induced in thorax and legs, 
particularly in flight muscle and leg muscle, whereas 
induction in the head and abdomen was either greatly 
reduced or absent. The GstD-deltaARE-GFP reporter 
was also observed to produce slightly more expression 
in young flies in the upper abdomen. The mean GFP 
intensity throughout the body was quantified from 
captured images of multiple flies using Image J 
software, and this analysis confirmed the up-regulation 
of both reporters during aging, in both males and 
females (Figure 2B). Despite the absence of GFP 
induction in head and abdomen tissue, the mean 
intensity of fluorescence produced by the GstD- 
deltaARE-GFP reporter throughout the fly was 
comparable to the un-mutated reporter (Figure 2B), as 
the expression in thorax and legs was relatively greater 
(see Figure 2A); and this difference may be due to some 
effect of the different chromosomal insertion sites on 
the overall expression levels for the reporters. Taken 
together, these data confirm the up-regulation of the 
GstDl gene during aging in the majority of adult 
tissues, and suggest that efficient expression in the head 
and abdomen may require the consensus ARE motif 
(Diagrammed in Supplemental Figure S4). 

DISCUSSION 

A core set of stress response genes shared with aging 

Changes in gene expression that were common to aging 
and all tested stresses identified a core set of stress 
response genes. Pepck encodes an enzyme critical in 
gluconeogenesis and glyceroneogenesis, and its up- 
regulation may be part of a basic metabolic adaptation 
to stress [19]. Interestingly, in mice, over-expression of 
PEPCK specifically in muscle tissue increases 
movement, life span and muscle mitochondrial 
proliferation [20]. Also among the core set of induced 
genes were Starvin, ref(2)P, and the Hsp genes Hsp70, 
Hsp83 and 1(2) efl. Starvin encodes a co-chaperone 
involved in autophagy and muscle maintenance and its 
up-regulation is consistent with its role in protein 
turnover and the cellular response to proteotoxicity [21]. 
Similarly, ref(2)P encodes a component of the 
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mitochondrial unfolded protein response (mUPR) 
pathway, consistent with proteotoxicity in the mito- 
chondrial compartment. Hsp genes are induced in 
response to protein denaturation and misfolding through 
activation of the HSF transcription factor, which in turn 
binds to HSEs in the Hsp gene promoters and activates 
transcription [1, 2]. In stressed cells certain Hsps have 
been shown to function to reduce proteotoxicity by 
favoring protein re-folding as well as the turnover of 
damaged proteins through the ubiquitin/proteasome and 
autophagy pathways [2, 22]. Induction of Drosophila 
Hsp genes during normal aging and upon oxidative 
stress has been shown to be dependent upon functional 
HSEs in the gene promoters, consistent with an 
increased abundance of misfolded proteins and conse- 
quent HSF activation under these conditions [4, 5]. The 



presence of Hsp genes in the set up-regulated by aging 
and each stressor indicates that protein 
denaturation/misfolding and HSF activation are 
common features of aging and each of the tested 
stresses. Induction of Hsps during aging may be part of 
a stress response that favors fly function by helping the 
fly to cope with aging-associated proteotoxicity [3]. 
Consistent with this idea, increased expression of 
certain Hsps is associated with increased fly life span 
[23-25]. In addition, it is also possible that chronic Hsp 
induction, particularly at late ages, may sometimes be 
mal-adaptive [2]. Interestingly, 8/18 of the common up- 
regulated genes and 21/32 of the common down- 
regulated genes have uncharacterized functions, 
indicating that there is much yet to be learned about the 
core stress response. 



Table 5. Features unique to aging 

a) GO enrichment terms for genes uniquely up-regulated in aging 



GO:0006952 


defense response(24) 


1.04E-10 


GO:0042742 


defense response to bacterium(17) 


2.52E-09 


GO:0009617 


response to bacterium(17) 


2.61E-08 


GO:0050830 


defense response to Gram-positive bacterium(l 1) 


1.35E-07 


GO:0051707 


response to other organism(19) 


1.29E-06 


GO:0009607 


response to biotic stimulus(19) 


1.55E-06 


GO:0006955 


immune response(18) 


4.36E-06 


GO:0002376 


immune system process(20) 


1.56E-05 


GO:0051704 


multi-organism process(21) 


7.84E-04 


GO:0019731 


antibacterial humoral response(8) 


0.002999 


GO:0009620 


response to fungus(7) 


0.010042 


GO:0006959 


humoral immune response(l 1) 


0.015651 


GO:0019730 


antimicrobial humoral response(lO) 


0.039404 


b) GO enrichment terms for genes uniquely down-regulated in aging 


GO:0006091 


generation of precursor metabolites and energy(53) 


6.12E-32 


GO:0045333 


cellular respiration(44) 


3.06E-29 


GO:0015980 


energy derivation by oxidation of organic compounds(44)) 


2.36E-27 


GO:0042773 


ATP synthesis coupled electron transport(29) 


9.71E-20 


GO:0006119 


oxidative phosphorylation(30) 


1.06E-19 


GO:0022900 


electron transport chain(30) 


1.09E-18 


GO:0022904 


respiratory electron transport chain(29) 


1.88E-18 


GO:0042775 


mitochondrial ATP synthesis coupled electron transport(27) 


4.70E-18 


GO:0055114 


oxidation-reduction process(62) 


8.74E-10 


GO:0006084 


acetyl-CoA metabolic process(17) 


1.56E-09 
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i~\ r\ . f\(\f\s AAA 

00:0006099 


tricarboxylic acid cycle(16) 


s~ cor? A A 

6.58b-09 


/"*/"V C\f\ A /~ ^ E / 

GO:0046356 


acetyl-CoA catabolic process(16) 


/ ror c\c\ 

6.58E-09 


CjU:0009109 


coenzyme catabolic process(16) 


1 AOr AO 

1.02b-08 


GO:0009060 


aerobic respiration(16) 


1.57E-08 


GO:0051187 


cofactor catabolic process(16) 


1.57E-08 


00:0006 /32 


coenzyme metabolic process(21) 


3.3 /b-0/ 


i~\ r\ . A A A/I 1 OA 

00:0006120 


mitochondrial electron transport, NADH to ubiquinone(13) 


^7 /^OT7 A^7 

7.68b-07 


rTk.AACI 1 O/C 

00:0051 loo 


cofactor metabolic process(22) 


1 11C A/C 

1.21b-06 


00:0016310 


Phosphorylation(47) 


1 O /I T7 A/T 

1.84b-06 


r*A.AAA^7m 

00:0006/93 


phosphorus metabolic process(53) 


/.4/b-06 


00:0006796 


phosphate-contammg compound metabolic process(53) 


7.47E-06 


r~\ aaa/^ooa 

GO:0006839 


mitochondrial transport! 13) 


O O 1 T -1 A A 

3.31E-04 


PA.AAA^I TO 

00:0006123 


mitochondrial electron transport, cytochrome c to oxygen(7) 


A AA1AAO 

0.002098 


rT\.AAAAAC^ 

00:0009056 


catabolic process(46) 


A AA1 1 1 1 

0.0031 13 


PA.AA/1 A TO 1 

00:0044281 


small molecule metabolic process(58) 


A A A 1 TOO 

0.003288 


AA-AAA^AAC 

00:0007005 


mitochondrion organization(15) 


A AAO C A A 

0.003544 


00:0006626 


protein targeting to mitochondrion(9) 


A AACAA/I 

0.005904 


AA.AA^ACOC 

00:0070585 


protein localization in mitochondrion(9) 


A AACAA/I 

0.005904 


PA.AATO^CC 

00:0072655 


establishment of protein localization in mitochondrion(9) 


A AACAA/I 

0.005904 


r*A.AAA'7TO') 

00:000/283 


Spermatogenesis(22) 


A A AT 1 /IT 

0.00/14/ 


00:0048232 


male gamete generation(22) 


A A AT O O /I 

0.007884 


GO:0043648 


dicarboxylic acid metabolic process(7) 


A A 1 A /*" O /*" 

0.010686 


GO:0006096 


Glycolysis(8) 


0.02534 


GO:0006122 


mitochondrial electron transport, ubiquinol to cytochrome c(6) 


0.026633 


GO:0044248 


cellular catabolic process(36) 


0.032559 


GO:0006006 


glucose metabolic process(lO) 


0.042632 



Each stress has unique gene expression changes 

While the pattern of gene expression changes during 
aging shared features with each of the tested stresses, 
each stress also had unique features (Summarized in 
Table 3). For example, hydrogen peroxide stress caused 
up-regulation of numerous genes involved in 
developmental pathways and signaling pathways, 
consistent with the fact that hydrogen peroxide also 
normally functions as a signaling molecule, during 
development and in adults, in Drosophila and other 
metazoans [26-30]. Genes uniquely up-regulated upon 
heat stress included ones of the Hsp60-class, which are 
important for protein trafficking and protein import into 
organelles [31]. Down-regulated genes unique to heat 
stress were enriched for the GO terms Defense response 
and Melanization defense response, suggesting that 
responses to wounding and bacterial challenge may be 
impaired. Finally, genes uniquely up-regulated in 



response to ionizing radiation included several 
proteasome subunit genes, which may indicate a 
particular requirement for protein turnover, perhaps in 
response to protein backbone cleavage, or alternatively 
this might reflect the critical role of the proteasome in 
DNA repair [32]. 

Aging has both shared and unique features relative 
to the tested stresses 

While aging shared features with each stress, aging was 
found to be more similar to the stresses most associated 
with oxidative stress (hyperoxia, hydrogen peroxide, 
ionizing radiation) than it was to heat stress. These 
observations are consistent with the conclusion that 
aging eukaryotic cells are in a pro-oxidant state [6, 33] 
associated with up-regulation of oxidative stress- 
response genes including ones encoding Gsts [34]. In 
addition to the shared features, a number of gene 
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expression changes were found to be unique to aging 
(Supplemental Table SI). For example, the gene 
encoding the mitochondrial form of superoxide 
dismutase {MnSOD or SOD2) was uniquely down- 
regulated during aging, and this is of potential interest 
given the fact that augmenting the expression of MnSOD 



can favor life span in adult flies [35, 36] and in C. 
elegans [37]. Also uniquely up-regulated during aging 
were the odorant receptor genes Obp56a and Obp57d, 
which is interesting in light of reports of negative 
effects of other odorant receptor genes on fly life span 
[38]. 



Table 6. Features common to aging and individual stress factors 

a) GO enrichment terms for genes up-regulated in aging and in hyperoxia 



GO:0009408 


response to heat(16) 


2.64E-11 


GO:0009266 


response to temperature stimulus(17) 


1.15E-10 


GO:0006950 


response to stress(34) 


9.95E-08 


GO:0035079 


polytene chromosome puffing(5) 


3.88E-05 


GO:0035080 


heat shock-mediated polytene chromosome( puffing(5) 


3.88E-05 


GO:0009628 


response to abiotic stimulus(17) 


8.38E-05 


GO:0044271 


cellular nitrogen compound biosynthetic process(14) 


0.00202 


GO:0034605 


cellular response to heat(6) 


0.002657 


GO:0019731 


antibacterial humoral response(7) 


0.009046 


GO:0009156 


ribonucleoside monophosphate biosynthetic process(5) 


0.032356 


GO:0009161 


ribonucleoside monophosphate metabolic process(5) 


0.032356 


GO:0006564 


L-serine biosynthetic process(3) 


0.04487 



b) GO enrichment terms for genes down-regulated in aging and in hyperoxia 



GO:0006508 


Proteolysis(56) 


6.01E-18 


GO:0045297 


post-mating behavior(7) 


9.33E-04 


GO:0008152 


metabolic process(131) 


0.00614 


c) GO enrichment terms for genes up-reguated in aging and in hydrogen peroxide 




GO:0009408 


response to heat(l 1) 


4.75E-06 




GO:0035079 


polytene chromosome puffmg(5) 


1.44E-05 




GO:0035080 


heat shock-mediated polytene chromosome puffmg(6) 


1.44E-05 




GO:0006950 


response to stress(27) 


2.86E-05 




GO:0009266 


response to temperature stimulus(l 1) 


7.87E-05 




GO:0034605 


cellular response to heat(6) 


8.32E-04 




GO:0044271 


cellular nitrogen compound biosynthetic process(13) 


0.001277 




GO:0009069 


serine family amino acid metabolic process(5) 


0.005345 




GO:0044281 


small molecule metabolic process(24) 


0.005367 




GO:0006564 


L-serine biosynthetic process(3) 


0.024923 




GO:0051707 


response to other organism(l 1) 


0.032029 




GO:0009607 


response to biotic stimulus(l 1) 


0.035404 





www.impactaging.com 



779 



AGING, November 2012, Vol.4 No.ll 



d) GO enrichment terms for genes down-reguated in aging and in hydrogen peroxide 



GO:0006091 


generation of precursor metabolites and energy(12) 


0.001493 


GO:0022900 


electron transport chain(9) 


0.001518 


GO:0055114 


oxidation-reduction process(22) 


0.002986 


GO:0045333 


cellular respiration(lO) 


0.003694 


GO:0015980 


energy derivation by oxidation of organic compounds(lO) 


0.008153 


GO:0022904 


respiratory electron transport chain(8) 


0.010626 


GO:0042775 


mitochondrial ATP synthesis coupled electron transport(7) 


0.043253 


e) GO enrichment terms for genes up-regulated in aging and in heat stress 


GO:0009408 


response to heat(18) 


1.04E-21 


GO:0009266 


response to temperature stimulus(18) 


1.68E-19 


GO:0009628 


response to abiotic stimulus(18) 


8.84E-13 


GO:0006950 


response to stress(22) 


1.70E-08 


GO:0006457 


protein folding(l 1) 


1.90E-07 


GO:0035079 


polytene chromosome puffmg(5) 


4.29E-07 


GO:0035080 


heat shock-mediated polytene chromosome puffing(5) 


4.29E-07 


GO:0034605 


cellular response to heat(6) 


1.27E-05 


GO:0001666 


response to hypoxia(6) 


0.001193 


GO:0070482 


response to oxygen levels(6) 


0.002157 


GO:0042221 


response to chemical stimulus(14) 


0.031137 


f) GO enrichment terms for genes down-regulated in aging and in heat stress 


GO:0006508 


Proteolysis(25) 


1.75E-07 


g) GO enrichment terms for genes up-regulated in aging and Ionizing radiation 


GO:0006950 


response to stress(33) 


2.17E-06 


GO:0009408 


response to heat(12) 


5.07E-06 


GO:0035079 


polytene chromosome puffmg(5) 


5.04E-05 


GO:0035080 


heat shock-mediated polytene chromosome puffmg(5) 


5.04E-05 


GO:0009266 


response to temperature stimulus(12) 


1.05E-04 


GO:0034605 


cellular response to heat(6) 


0.003611 


GO:0009069 


serine family amino acid metabolic process(5) 


0.018171 


GO:0044271 


cellular nitrogen compound biosynthetic process(13) 


0.022327 


GO:0033554 


cellular response to stress(18) 


0.024714 
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h) GO enrichment terms for genes down-regulated in aging and Ionizing radiation 



GO:0006091 


generation of precursor metabolites and energy(17) 


2.30E-07 


GO:0045333 


cellular respiration(14) 


1.95E-06 


GO:0015980 


energy derivation by oxidation of organic compounds(14) 


6.16E-06 


GO:0022900 


electron transport chain(l 1) 


4.47E-05 


GO:0055114 


oxidation-reduction process(26) 


2.65E-04 


GO:0006119 


oxidative phosphorylation(lO) 


2.91E-04 


GO:0022904 


respiratory electron transport chain(lO) 


2.91E-04 


GO:0042775 


mitochondrial ATP synthesis coupled electron transport(9) 


9.52E-04 


GO:0042773 


ATP synthesis coupled electron transport(9) 


0.00164 



The aging gene expression pattern indicates a failure 
in mitochondrial maintenance 

The changes in gene expression that were found to be 
unique to aging included up-regulation of numerous 
innate immune response genes, and down-regulation of 
numerous mitochondrial metabolism genes, including 
ones encoding components of the ETC (Table 4). 
While up-regulation of innate immune response genes is 
a feature of aging that is shared with hyperoxia [6] 
(Supplemental Table S3), the number of up-regulated 
innate immune response genes was significantly greater 
for aging, resulting in many changes in this category 
that were unique to aging. Similarly, down-regulation 
of mitochondrial and ETC genes is a feature of aging 
that is shared with ionizing radiation and hydrogen 
peroxide stress (Supplemental Table S3), but the 
number of down-regulated mitochondrial and ETC 
genes was greater for aging, resulting in many changes 
in this category that were unique to aging. Girardot et 
al [39] examined gene expression changes during 
Drosophila aging separately for the head, thorax and 
abdomen, and found that down-regulation of 
mitochondrial genes is observed preferentially in the 
thorax; because the thorax is composed primarily of 
flight muscle this observation suggests that 
mitochondrial gene down-regulation may occur 
preferentially in muscle tissue. 

Up-regulation of innate immune response genes during 
Drosophila aging is in part due to a dramatic increase in 
microbial load during aging, as eliminating bacteria 
reduces the response [13]. However, innate immune 
response genes are still up-regulated during aging in the 
absence of detectable microbes, suggesting additional 
mechanisms for activation of these genes during aging. 
Consistent with this conclusion, innate immune 
response genes are also up-regulated in response to 



oxidative stress caused by hyperoxia ([6]; this study), 
and therefore one possibility is that an aging-related 
failure in mitochondrial maintenance leads to oxidative 
stress that can induce innate immune response gene 
expression. Similarly, studies in mammals reveal that 
damaged mitochondria also release DNA fragments and 
formyl-peptides that can induce innate immune 
response genes [40], and therefore this may be an 
additional mechanism for innate immune response gene 
induction during aging that is a consequence of a failure 
in mitochondrial maintenance. The across-the-board 
down-regulation of Drosophila mitochondrial genes, 
ETC genes and mitochondrial metabolism genes 
observed during aging suggests a possible mechanism 
for a failure in mitochondrial maintenance during aging 
(Diagrammed in Supplemental Figure S5). The ETC 
and mitochondria turn over at a basal rate, and more 
rapidly in response to signals such as starvation, and a 
reduced rate of replacement is expected to result in 
longer-lived structures that will be more susceptible to 
time-dependent damage and malfunction. This idea is 
consistent with the observed accumulation of structural- 
ly abnormal mitochondria during Drosophila aging [41- 
44], reduced mitochondrial transcription [45], decreased 
ATP and increased production of ROS [46]. Decreased 
ATP flux is expected to reduce rates of bulk protein 
synthesis and turnover, and increased ROS will increase 
protein damage, consistent with the accumulation of 
damaged and misfolded proteins (proteotoxicity) and 
the induction of Hsp genes [2, 22, 24]. 

Taken together, the data support a model wherein the 
down-regulation of mitochondrial and ETC genes 
during aging leads to a failure in mitochondrial 
maintenance and the accumulation of abnormal 
mitochondria, which in turns leads to oxidative stress 
and proteotoxicity; these stresses in turn cause the 
oxidative-stress-like and proteotoxic-stress-like patterns 
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of gene expression observed during aging 
(Supplemental Figure S5). Placing oxidative stress 
down-stream of an aging-associated failure in 
mitochondrial maintenance is consistent with the 
observation that oxidative stress correlates with, but 
does not directly regulate life span in Drosophila [47], 
and with the implication of mitochondrial malfunction 
in mammalian aging-related metabolic disorders [48]. 
Consistent with the importance of mitochondrial 
maintenance in aging, certain interventions that increase 
mitochondrial proliferation, such as over-expression of 
PGC1 alpha in gut tissue, have recently been reported to 
increase life span and tissue function in aging 
Drosophila [49, 50], and PGC1 alpha activity is also 
implicated in maintaining tissue function during aging 
in mammals [51]. In contrast, other manipulations that 
increase Drosophila mitochondrial proliferation, such as 
increased tissue-general expression of PGC1 alpha [50] 
or cyclin D/Cdk4 [52] had negative consequences for 
life span and oxidative stress levels, and taken together 
these studies indicate that effective interventions in 
mitochondrial maintenance during aging will require 
tissue-specific targeting. Notably, certain carefully- 
timed interventions that reduce activity of ETC 
components have been shown to increase life span in 
both invertebrates and mammals [53-55], and this might 
function through a hormetic response to increase 
production of new mitochondria, or conceivably by 
inhibiting the activity of abnormal mitochondria. 
Critical questions for the future include determining the 
causes and mechanisms for the observed down- 
regulation of mitochondrial and ETC genes during 
aging - a pattern shared by Drosophila and mammalian 
tissues [6, 9]. Possible explanations include the 
inherently shorter-lived nature of mitochondrial genome 
sequences relative to nuclear genome sequences, genetic 
conflicts resulting from the uni-parental inheritance of 
mitochondrial genomes, and trade-offs between the 
costly production of new mitochondria and investments 
in growth, sexual differentiation and reproduction [56- 
61](Supplemental Figure S5), and these will be 
interesting areas for future research. 

METHODS 

Drosophila culture, microscopy and stress treatments. 
Drosophila melanogaster flies were cultured on a 
standard agar/dextrose/corn meal/yeast media at 25°C 
[62]. The transgenic strains GstDl-GFP and GstDl- 
deltaARE-GFP were generously provided by Dirk 
Bohmann [15]. Age-synchronized cohorts of flies were 
generated by collecting newly-eclosed flies over a 
period of 48 hours, followed by maintenance in vials at 
approximately 20 flies per vial, with every-other day 
transfer to fresh media, until the indicated age time 



points. Visible images, GFP fluorescence images, and 
image overlays for flies were generated using the Leica 
MZFLIII fluorescence stereomicroscope. GFP fluore- 
scence was quantified using captured GFP images and 
Image J software, with mean and standard deviation 
calculated using 6 flies per sample. Flies used for 
stress treatments, RNA analyses and microarray 
analyses were generated as follows: males of wild-type 
strain Oregon-R were crossed to virgins of transgenic 
laboratory stock w[1118] ;rtTA(3)E2/TM3 Sb to 
generate hybrid progeny of genotype 
w[1118];rtTA(3)E2/+, as was used for previous 
microarray analyses [63], and 9-10 day-old male flies 
were used for each stress treatment. Old flies were 61 
days of age, which corresponds to approximately the 
50% survival point for the cohort [6]. Vials containing 
1% sucrose were prepared by adding 1.5 ml of 1% 
sucrose in deionized water to a Drosophila vial 
containing a single folded Kimwipe (Kimberly-Clark). 
For each stress treatment and the sugar-treated controls, 
replicate vials of 25 flies each were subjected to the 
treatment, and then the flies from each vial were 
separately processed for RNA, and each sample was 
used to generate probe for one micro-array 
hybridization, such that each treatment is represented by 
at least three biological replicates. For hyperoxia 
treatment flies in standard food vials were subjected to 
100% oxygen atmosphere for 5 days as previously 
described [6]. For ionizing radiation treatment flies in 
standard food vials were irradiated with 5666Rads/hour 
for 16 hours using a Cesium source (Grammacell 40- 
Cesium 137, Atomic Energy, Canada) at the USC 
Norris Cancer Center facility, and then transferred to 
1%> sucrose vials for two days followed by processing 
for RNA. For qPCR analysis 9 hour irradiation samples 
were also generated. Because ionizing radiation is 
inhibitory to transcription, the two-day recovery period 
was included to allow the gene expression response to 
develop; recovery in sucrose vials was employed 
because the newly-irradiated flies have greatly reduced 
mobility and will adhere to the surface of a regular food 
vial. For hydrogen peroxide treatment flies were placed 
in sucrose vials adjusted to 3% hydrogen peroxide for 
two days, and then processed for RNA. For heat stress 
treatment flies were placed in sucrose vials at 37°C for 
5.5 hours and then processed for RNA. Controls for the 
effects of sucrose vials ("sugar-treated controls") were 
generated by placing flies in sucrose vials for two days 
prior to processing for RNA. 

RNA isolation and microarray hybridization. An 
average of 35 jug RNA was isolated from groups of 25 
adult male flies using Trizol reagent (Life Technologies, 
Grand Island) according to the manufacturer's 
instructions. The RNA was further purified using the 
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RNAqueous kit, and concentration was determined 
using NanoDrop spectrophotometer. A portion of the 
RNA (3 jug) was fractionated on 1.0% agarose gels to 
determine purity. 10 jug of total RNA was then used as 
substrate to generate biotinylated cRNA according to 
standard Affymetrix protocol (Childrens Hospital, Los 
Angeles, CA)[6]. A total of 35 Affymetrix gene chips 
were analyzed including at least four biological 
replicates for each experimental condition and control, 
with the exception of heat stress in which one array was 
omitted due to poor quality. The old, hyperoxia, and 
young samples were derived from our previous study 
[6] in which six arrays were used for the hyperoxia and 
young conditions and four arrays were used for old 
flies. Quantitative real-time RT-PCR analyses [64] and 
Northern blot analyses [65] were performed as 
previously described, using RNA samples derived 
independently from those used for the microarrays. 

Statistical analysis of microarray data. Gene expres- 
sion measures were computed based on a non-linear 
multi-chip model of the perfect match signal [66]. This 
approach enables the separation of specific and non- 
specific components of the microarray signal and 
circumvents the issue of saturation bias in the high- 
intensity range. The background and concentration 
parameters were both fit within a single global routine 
(rather than estimating the background parameter before 
computing gene expression measures), and the model 
that best described the observed data selected. Linear 
modeling and empirical Bayes analysis [67] was 
performed in the R statistical programming language 
( http ://www.r-proj ect.org/ ) using the Limma: Linear 
Models for Microarray Data package [67] to identify 
genes significantly differentially expressed in response 
to multiple stressors or during aging; Limma computes 
an empirical Bayes adjustment for the t-test. Because 
the identification of genes altered in multiple conditions 
was a major objective of this study, a nested F-test 
approach was employed as this can be more powerful at 
detecting genes altered in multiple contrasts. Multiple 
testing was corrected for using the Benjamini and 
Hochberg method, which controls the false discovery 
rate (FDR) [68] in this framework on a per-gene basis 
(but not across contrasts). Using this robust method, 
genes were found to be significantly differentially 
expressed both by biological and statistical criteria 
(±1.2 fold change, FDR 1% (p < 0.01); Supplemental 
Table S4); Gene expression changes of ±1.5 fold were 
used for subsequent comparisons, as indicated. 
Hierarchical cluster analysis of the top 1000 
differentially expressed genes for each condition based 
on the F-test /?-value from the linear model fit was 
performed to visualize the gene expression patterns 
across different stressors, using the R package mclust. 



The microarray data discussed in this study have been 
deposited in the National Center for Biotechnology 
Information Gene Expression Omnibus (GEO) [69] and 
are accessible through GEO Series number GSEXXX. 

Functional annotation and statistical overrepresentation 
of Gene Ontology classifications. Statistically over- 
represented GO categories were identified using 
Flymine [70], by the calculation of a p-value denoting 
the probability that the observed numbers of counts 
could have resulted from randomly distributing a 
particular GO term between the test and the reference 
group. Multiple testing was controlled for using the 
Holm-Bonferoni method. 

Statistical significance of overlapping gene sets. The 
statistical significance of the overlap between various 
gene sets was evaluated by computing the ^-value 
representing the probability of obtaining the observed 
number of overlaps by chance under a hypergeometric 
distribution, using the R function phyper [71]. 

Identification of enriched GO terms and corrections for 
effects of sucrose vials. Gene annotations for the 
AffyDrosGenomel arrays were updated to the latest 
information from Flybase using the online tool Flymine 
[70] for all genes with expression altered >1.5 fold. 
Gene annotations identified by Flymine as matching 
more than one entry in the current database were 
resolved where possible, as follows: The Affymetrix 
probe ID was obtained from the limma files, and the 
corresponding probe sequence was obtained from the 
Affymetrix website. The probe sequence was then used 
to query the current Drosophila genome sequence 
annotation using the Flybase website and BLAST 
function to identify the correct gene. Ten probe 
sequences had ambiguous match that could not be 
resolved and were not included in the GO term analyses 
(Supplemental Table S5), and an additional 27 
identifiers did not match genes in the current database. 
To control for any possible effects on gene expression 
patterns caused by two days maintenance of flies in 
sucrose vials, GeneChip analysis was performed on flies 
transferred to sucrose vials for two days in the absence 
of added stressors as a control. 258 genes were found to 
be up-regulated and 362 genes were found to be down- 
regulated relative to controls maintained on normal 
media (Supplemental Table S6), and these gene sets had 
no GO terms enriched among the up-regulated genes, 
and 4 GO terms enriched among the down-regulated 
genes: Proteolysis, Post-mating behavior, Insemination, 
and Lipid metabolic process (Supplemental Table S7). 
The genes that were up-regulated and down-regulated in 
response to sucrose were subtracted from the lists of 
genes up-regulated and down-regulated in response to 
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hydrogen peroxide and ionizing radiation treatment to 
generate the final lists presented in Tables 1, 3, 4, and 
Supplemental Tables SI -S3. While this simple 
subtraction procedure does not account for possible 
gene expression changes caused by interactions of 
sucrose with the stressors, we observe that the major 
GO term categories enriched in the gene sets up- 
regulated and down-regulated by hydrogen peroxide 
and ionizing radiation do not differ significantly when 
the gene expression changes caused by sucrose alone 
are included or excluded from the analysis (compare 
Supplemental Tables S8 and S9 where the effects of 
sucrose are included, to Supplemental Tables S10 and 
Sll where the effects of sucrose are excluded). In 
addition, cluster analysis demonstrated that the gene 
expression changes caused by the hydrogen peroxide 
stress treatment and ionizing radiation stress treatment 
were more similar to each other and to aging and 
hyperoxia than they were to the sucrose-treated control 
flies (Supplemental Figure S3), providing further 
evidence that the gene expression changes due to 
sucrose transfer do not make a significant contribution 
to the gene expression changes observed in the 
hydrogen peroxide and ionizing radiation samples. 
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Figure SI. Northern analysis of gene expression in response to selected stresses. Total RNA was isolated from 
flies subjected to each of the indicated stresses, and 4ug (IX) and 8ug (2X) quantities were analyzed for each sample, as 
indicated. The blot was hybridized successively with the indicated gene-specific probes, and visualized by autoradiography. 
Signals were quantified using the phosphorimager, and the fold induction relative to control was calculated for both the IX 
and 2X samples, the two values were averaged, and the result is presented below the lanes. For the Hsp70 gene-specific 
probe, both a darker exposure and lighter exposure are presented, as indicated. H202, hydrogen peroxide. 
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Color Key 




Figure S2. Cluster analysis and heatmap. A heatmap visualization of the gene expression levels across each condition for 
the top 1000 differentially expressed genes (according to the F-test p-value from the linear model fit). XRay, ionizing radiation. 
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Figure S3. Gene RNA levels in response to selected stresses. Quantitative real-time RT-PCR analysis was used to 
determine RNA levels for the genes Hsp22, Drosomycin, Ade3 and CG11089 in response to selected stresses, in both male 
and female flies, as indicated. HS, heat stress; IR, ionizing radiation. Stress treatment RNA levels were compared to 
control using unpaired, two-sided t-tests, and statistically significant differences (p < 0.05) are indicated with asterisk. 



www.impactaging.com 



788 



AGING, November 2012, Vol.4 No.ll 



- GstD2 



Foxo ARE 



TRE 



+1 



GstDl — 



GstD-GFP 



ttgtttac 



Foxo 



TCAgcATGACcggGCAaaaa 



ARE 



1241 


-869 




tatgagtca 









.+1 



GFP - 



TRE 



GstD-deltaARE-GFP 



r> +1 





ttgtttac 




TCAgcACTCGaggGCAaaaa 




tatgagtca 









GFP 



Foxo 



deltaARE 



TRE 



Figure S4. Diagram of control and mutant GstDl-GFP reporters. Diagram of the genomic region containing 
the GstD2 and GstDl genes and the structure of the GstD-GFP and GstD-deltaARE-GFP reporter constructs. The 
conserved Foxo and ARE motifs are indicated; red lettering indicates the bases mutated in the deltaARE construct. 
Numbering is relative to the +1 start site for transcription of the GstDl gene; the diagram is not to scale. 
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Figure S5. Model for gene expression changes during aging. The 

genes that comprise the Insulin/IGFl-like signaling pathway (IIS), the TOR 
signaling pathway (TOR) and the sex-determination pathway (SD) exhibit 
antagonistic pleiotropy (AP) and sexual antagonistic pleiotropy (SAP), in 
that these pathways act to promote growth, differentiation, sexual 
differentiation and reproduction at the expense of mitochondrial (Mito) 
gene expression and mitochondrial turnover. Decreased mitochondrial 
turnover results in the accumulation of abnormal and malfunctioning 
mitochondria, and in turn decreased ATP flux and increased ROS 
production. Decreased protein turnover due to decreased ATP flux and 
increased protein damage due to increased ROS production combine to 
increase levels of protein damage and misfolding, leading to 
proteotoxicity. Release of mitochondrial DNA fragments, formyl peptides 
and ROS from abnormal mitochondria (combined with increased microbial 
load) cause dramatic induction of innate immune response genes, 
including anti-microbial peptide (AMP) genes. The increased ROS 
production induces additional types of oxidative stress response genes, 
including Gluthione-S-transferase (Gst) genes. Protein misfolding activates 
the heat shock transcription factor (HSF) and mitochondrial unfolded 
protein response (mUPR) to induce proteotoxicity-response genes, 
including heat shock protein (Hsp) genes. Classes of genes up-regulated 
during aging are indicated with red color font. 
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